Crumpling a Thin Sheet 
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Crumpled sheets have a surprisingly large resistance to further compression. We have studied the 
crumpling of thin sheets of Mylar under different loading conditions. When placed under a fixed 
compressive force, the size of a crumpled material decreases logarithmically in time for periods up 
to three weeks. We also find hysteretic behavior when measuring the compression as a function 
of applied force. By using a pre-treating protocol, we control this hysteresis and find reproducible 
scaling behavior for the size of the crumpled material as a function of the applied force. 

PACS numbers: 68.60.Bs, 89.75.Da, 62.20.Fe, 46.32,+x 



If you, the reader, were to rip a page from this jour- 
nal and crumple it, squeezing it with your hands into a 
ball as hard as you can, the resulting object is still more 
than 75% air. What gives this crumpled sheet its sur- 
prising strength and how does the ultimate size of the 
sheet depend on the forces applied? 

The energy stored in crumpled sheets has been inves- 
tigated in the limit where the thickness, 8, of the sheet 
is much less than the lateral dimension, L (e.g., the di- 
ameter of a circular sheet) |[ ||, [§ Hi- Analysis of 
these objects has focused on two characteristic structures 
caused by the crumpling: singular, conical points and the 
curved ridges, which store most of the energy, connecting 
them. As the sheet is crumpled farther, the ridges must 
collapse so as to remain within the confining container; 
and an increasing number of smaller ridges must form. 
In the limit where the diameter of the container, D, is 
much smaller than the initial size of the sheet but still 
much larger than the sheet thickness (8 <C D <C L) a 
scaling relation is predicted between the energy stored 
and the size of the resulting crumple: E oc D@ . This pre- 
diction assumes that the forces resisting crumpling are 
conservative. One knows however, that frictional forces 
are involved as the sheet rubs against itself and against 
the constraining walls. How do these forces affect the re- 
sulting energy balance? In addition plastic flow can occur 
in regions of high curvature effectively constraining the 
geometry of the energy bearing structures. 

The notion of scaling properties describing a crushed 
sheet goes back to the initial experiments of Kantor et 
al. []5| Further investigations of such sheets were made 
by Gomes and collaborators || Q . More broadly, scaling 
behavior of force with compression in tenuous structures 
was reported for random fibrous material by Baudequin 
et al. H and for colloidal aggregates by several groups || . 
Here we report distinctive forms of scaling not seen in 
those previous studies. 

To study how the external dimensions of a crumpled 
sheet depend on the confining force, we have placed a 
large circular (L = 34cm) sheet of thin (8 — 12.5/xm) 
aluminized Mylar under a weighted piston in a plas- 
tic cylindrical cell of diameter D = 10.2cm. The cylinder 



is surrounded by a box to protect it from air currents 
in the room and is mounted on an optical table to mini- 
mize the effects of vibration. We measure compression by 
measuring the height, h, of the piston above the base. No 
special care is taken in the initial crumpling that allows 
the large sheet to be placed within the cylinder. 

After a mass, M, is placed on the piston, the piston 
settles to a new height held up by the crumpled sheet 
below it. However, this behavior is not as simple as one 
might have expected: h continues to decrease long after 
the mass is initially introduced. As shown in Fig. 0(a), 
over seven decades in time, 2 x 10 _1 sec <t<2x 10 6 sec, 
the entire duration of the experiment, the Mylar height 
decreases logarithmically: 



h = a — b\og(t/ sec) 



(1) 



where a and b are constants. Even after three weeks the 
piston did not reach its asymptotic height. Similar relax- 
ation occurred in other tenuous materials, as shown in the 
insets to Fig. 0(a). Absorbent tissue paper |ll[] showed a 
logt dependence like that of the Mylar. Cotton balls 
showed a slightly smaller relaxation over several weeks, 
in which the logarithmic dependence extended only up to 
10 4 sec. Although this logarithmic behavior is robust and 
is found each time the Mylar is compressed, the constants 
a and b fluctuate between runs with no monotonic evo- 
lution of those constants as the sheet is crumpled more 
times. Fig. 0(b) shows a weak correlation between a and 
b for the same 200g load on the piston. 

Albuquerque and Gomes studied stress relaxation 
of crumpled aluminum foil under fixed strain and found 
stretched-exponential behavior quite different from the 
logarithmic dependence of strain that we measured at 
fixed stress. Aluminum foil is much more maleable than 
the Mylar we have used and presumably has more plastic 
flow along the ridges and vertices. This may be the origin 
of the different relaxation kinetics. Menon jl3| has found 
similar log t behavior as our results in stress relaxation. 

Relaxation continues even if the compression is 
stopped. Fig. ||(a) shows the time dependence of the 
height of a Mylar sheet after a 200g mass has been place 
on it. After 500sec. the piston is lifted very slightly and 
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FIG. 1: (a) Time-dependence of the height of a Mylar sheet 
compressed by a mass M = 200<?. The straight line is a fit to 
Eq. hi with a — 8.4cm and b = 0.14cm. Insets show the time 
dependence of the height for porous tissue and for cotton balls 
under the same conditions: M = 200(7. (b) The constants a 
and b from Eq. |l| for a series of runs under nominally identical 
conditions. No obvious evolution of the values were found as 
the sheet was crumpled successively more times. 



fixed in place so that it can no longer compress the ma- 
terial. After lOOOsec, the piston is released and quickly 
sinks to a new position comparable to where it would 
have been had it never been restrained; thereafter it re- 
laxes at nearly the same logarithmic rate as it did ini- 
tially. 

What is responsible for the slow decrease in height over 
these extended periods? External vibrations affect the 
compression. Fig. ||(b) shows h(t) both before and after 
the container is vibrated laterally with an acceleration 
of 2.5m/ sec 2 at 7.5Hz. Initially, the container was at 
rest and the logarithmic slope, 6, was determined. Af- 
ter approximately 200sec, vibrations were begun and b 
nearly doubled. This change, while substantial, is not so 
large as to indicate that the slope in the quiescent state 
is caused by the ambient vibrations in the room which 
are less than 10~ 2 m/sec 2 . 

This logarithmic decay of h suggests activated relax- 
ation (with an effective temperature T e g introduced to 
account for external vibrations JI3]) of the sheet between 
metastable minima separated by energy barriers, AE: 
dh/dt = exp[—AE / (kBT e ff)]. As h decreases, the 
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FIG. 2: (a) Time-dependence of the height of a compressed 
Mylar sheet before and after the compression was stopped. 
The initial time-dependence of the height of a Mylar sheet 
after a 200g mass has been place on it is shown by the filled 
squares. After approximately 500sec. the piston is lifted and 
fixed in place. After lOOOsec, the piston is released and the 
subsequent relaxation is shown by the filled circles, (b) h(t) 
for a compressed Mylar sheet before and after external vibra- 
tions are introduced. Initially the container was at rest. After 
approximately 200sec, lateral vibrations with an acceleration 
of 2.5m/ sec 2 at 7.5Hz were begun. The resulting logarithmic 
slope, b, nearly doubled. 



crumpled sheet becomes progressively stiffer suggesting 
that the barriers are increasing. To lowest order (i.e., 
AE w E - Ch) we find: 



dh/dt = -Cex P [-(£ - Ch)/(k B T eff )} 



(2) 



which has a solution of the form of Eq. [j] where a = 
di+Eo/C-iksTeff/C) HCX/k B T eff ) and b = k B T eff /C 
where d\ is a constant. As T e g increases (e.g., due to ex- 
ternal vibrations) the logarithmic slope, b, also increases 
as observed in the experiments. If the precise initial 
crumpling determines Eq and C, then a and b should 
be correlated between runs: a = d± + b(lnb + d 2 ) where 
d 2 = E /kBT e ff — InC- The curve in Fig. |l|(b) shows the 
fit with di = 6.2 and d^ — 9.7. We note however, that d\ 
and g?2 may also vary between trials. 

So far, we have reported the time dependence of the 
height, h, under a fixed weight. We are also interested 
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in determining the effect of the piston mass, M, on the 
height. There are two problems which interfere with ob- 
taining such data in a straight-forward manner. First, as 
we have emphasized, h varies logarithmically with time 
and it is thus impractical to find the final height of a 
crumpled sheet under a given load. Second, there is large 
hysteresis so that different loading procedures produce 
different values of h. Such hysteresis is already suggested 
by the variation between runs of the constants, a and b, 
describing the relaxation. Despite these problems we can 
obtain reproducible results by following two protocols for 
handling the crumpled sheets. 

In order to avoid the problems created by the slow 
relaxation, we measure the height after a fixed interval 
of time (e.g., lOOsec) and use this value, instead of the 
unattainable asymptotic (t = oo) value, for comparison 
between different loadings. Such a procedure means that 
we neglect relaxation processes slower than lOOsec. To 
obtain meaningful results we confine ourselves to com- 
paring the piston height for relatively large changes in 
the mass, AM so that AHam (the change in height mea- 
sured at t — lOOsec due to the increased piston mass) 
is much greater than Ah^t (the change in height due to 
waiting extra time after the mass has been added, e.g., 
between lOOsec and lday ): Ah am 3> Ah At- In what 
follows we measure hioo(M), the height at t = lOOsec 
after a mass, M, has been added to the piston. 

Dealing with the second problem of hysteresis requires 
a protocol for preparing a sheet prior to measurement. 
Hysteresis is evident if we start with a small initial mass 
and measure (as prescribed above) the height after sub- 
sequent increments AM have been added until a large 
final mass is applied. If we then remove this load and re- 
peat the measurement with the initial mass, we find that 
the newly measured height is smaller than that initially 
measured for the same mass. 

In order to avoid this hysteresis, we have followed the 
protocol of Baudequin et al. for studying compressed 
glass wool [||. We pre-train a sheet by first crumpling 
it under a mass Mi for a period of one day. The mass 
is then removed and we measure /iiqo(M) for M < Mi 
(as shown for Mi = 2.6kg by the open O's in Fig. |](a)). 
This curve is reproducible (shown by the X's) as long 
as Mi is not exceeded. However, if we exceed Mi and 
place a larger mass on the piston, we no longer obtain the 
same curve. With this larger mass M2 = 4.5fc<7 we obtain 
a new curve (shown by the solid triangles) that is itself 
reproducible as long as M 2 is not exceeded. We continue 
this measurement to obtain a family of curves each of 
which corresponds to a pre-training mass Mj. For each 
curve, considered separately, we find a scaling relation: 



h 100 (M,Mi) 



cM' 



(3) 



where hi is the asymptotic value of the height for the 
pre-treating mass Mj. Such a family of curves is shown 
in Fig. |(b). We find a = 0.53 ± 0.04 for all pre-treating 
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FIG. 3: (a) Hysteresis in the height versus mass. A crumpled 
Mylar sheet is first placed under a mass Mi = 2.6kg for one 
day. The mass is then removed and /iioo(M) is measured 
for M < Mi (open O's). If Mi is not exceeded, this curve 
is reproducible (X's). If Mi is exceeded, a different curve 
hioo(M) is found (solid triangles), (b) log(ftioo(M, Mi) — hi) 
versus logM for five separate pre-treating masses Mi. Solid 
lines are fits to Eq. ^ with a = 0.53. Dashed line shows the 
predicted h(t). Inset shows scaling between the asymptotic 
height, hi, and the pre-treating mass, Mi. 



masses. This is smaller than the value a = 0.67 found by 
Baudequin et al. || for fibrous glass wool. The inset to 
Fig. |3|(b) shows a scaling relation between the asymptotic 
height, hi, and Mf. hi oc M~ 7 with 7 w 0.8. Despite 
the problems associated with the logarithmic relaxation 
and hysteresis, Fig. ||(b) shows that it is possible, under 
well-defined conditions, to obtain reproducible data on 
crumpled sheets. 

We can predict the scaling relation between M and h if 
we assume a small volume fraction and a uniform network 
of ridges between nearly flat facets of characteristic size 
X that depends on compression. The energy of each 
facet is proportional to n{X/ S) 1 ^ 3 , where k is the bending 
modulus of the sheet and S is its thickness [| ||, |^, 
M. Facets are oriented at random, so there is roughly 
one facet per volume X 3 . Thus the energy contained 
in the volume V is E oc V/ X 3 k(X / 5) 1 / 3 . The space 
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occupied by one facet is large enough to contain a stack 
of X/S facets so that = 5/X is a volume fraction. 
Macroscopically, cf> = (n5(L /2) 2 ) / '{nh{D /2) 2 ), where h 
is the height of the container, D is its diameter and L is 
the diameter of the sheet. E can be re-expressed: 

E = kV/S 3 (X/S)-^ 3 = (K/6 3 ){nh(D/2) 2 )(j) 8 / 3 oc /T 5/3 . 

(4) 

On the other hand, the energy E supplied by the applied 
force Mg is given by E = Q M(h)gdh w Mgh. As noted 
previously, we have assumed here that all the work goes 
into ridge energy. We conclude that Mgh w h~ 5 ^ 3 : so 
that h oc M~ 3 / 8 . This predicted value of a — 0.375 is to 
be compared with the measured a = 0.54. A more careful 
estimate |H| yields h ks 0.31cm when M = 0.1kg. The 
predicted height, shown by the dashed line of Fig. ||(b), 
is a factor 1.5 to 3 lower than those measured for the 
smallest masses, and approaches the measured heights 
for the largest masses. 

The computer simulations of Lobkovsky et al. |Q in- 
dicated that the asymptotic scaling worked well when 
the volume fraction was no more that a few tenths of a 
percent. However, the experiments are well outside this 
regime, where the volume fraction is large (over one half). 
The size of the observed facets is far from uniform, and 
this nonuniformity could alter the stored energy. Fur- 
ther, we neglect friction which is expected to strengthen 
the crumpled sheet by a factor of order unity. Plastic 
deformation is expected to weaken it by at most a fac- 
tor of order unity. Considering these differences between 
the experiment and the theoretical model, the agreement 
between the two exponents is surprisingly good. 

The data presented here suggests that the answer to 
our question in the first paragraph, about how the ulti- 
mate size of a sheet depends on the forces applied, has a 
complicated answer. The logt dependence found in the 
relaxation suggests that there is dissipation, either due 
to friction or to plastic flow in regions of large curvature. 
It is unclear to what extent an analysis based only on 
conservative forces, which does not consider such effects, 
can fully describe the energetics of crumpling in real ma- 
terials. Despite that caveat, our data nevertheless do 
support a scaling picture of the crumpling process: after 
a sheet has been appropriately pre-treated, the height 
approaches an asymptotic value as a power-law. 

The results suggest several extensions such as study- 
ing how the exponent a depends on size and thickness in 
different materials. The most pressing question is the ori- 
gin of the logarithmic relaxation. We do not know if the 



behavior occurs because of plastic flow in the material 
or because of frictional sliding at the contacts between 
sheets. Neither effect is included in the theoretic or sim- 
ulation studies of crumpling. In order to determine the 
role of plastic flow in regions of large curvature, one could 
measure the time-dependent relaxation in rubber sheets 
where such plasticity would be minimized. 
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in preparing this manuscript. 
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